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Abstract 

Starting with the S-wave radial equation for an attractive central 
potential V(r), we give results for the n(principal quantum number) 
and the [i (reduced mass) dependence of i?„o(0), the S-wave radial 
waveunction at the origin, for potentials with definite curvature. 

1 Introduction 

Discovery of quark-antiquark atoms like charmonium in 1974 led to general 
investigations [1 J of the Schrodinger equation with a central potential V(r) 
representing the gq-potential. The motivation was to obtain results based 
on general properties of V{r) like its shape, since its precise form was then 
(and still is) unknown. 

Some results were obtained for the S-wave (£ = 0) bound state radial 
wave function R n o(r) [2]. Specifically, it was shown that i?2o(0) is larger 
(smaller) than i?io(0) provided V(r) was everywhere convex, V"(r) > 
(concave, V"(r) < 0). Prime denotes derivative with respect to r. The 
variation of R\q(0) with the reduced mass \i was also related to the curvature 
of the potential. This was directly proved from the radial equation for n = 1 
[3]. Both types of results mentioned above were proved for large n using the 
WKB approximation [3]. 

In this note we show that both types of results follow directly from the 
S-wave radial Schrodinger equation for all n. For notational simplicity, define 

S n (0) = [Rn0(0)} 2 . (1) 
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We will prove that for n = 1, 2, . . . , 

(a) if V'(r) > and V"(r) = for all r then 



S n (0) - S n+1 {0) 



d 

and — — 



~S n (0) 
A 4 



(b) if V'(r) > and V"(r) < for all r and K'(oo) is finite then 



5„(0) -5 n+ i(0) > and |- 



-5 n (0) 
A 4 



>0; 



(c) if V'(r) > and V"(r) > for all r and V'(0) is finite then 



5„(0)-5 n +i(0) <0 and 



5/x 



/i 



-5 n (0) 



< 0. 



(2) 



(3) 



(4) 



Case (a) corresponds to an attractive linear potential. This case is exactly 
solvable. Indeed, there are well known exactly solvable examples for the 
concave (Coulomb potential) and convex (harmonic oscillator) cases which 
satisfy the above inequalities. Explicit solutions of convex power law po- 
tentials r k with k > 1 and the concave log(r) potential satisfy the above 
inequalities [5]. With all this evidence at hand we believe that the above 
inequalities are really theorems. In the next section we present a result for 
non-zero I followed by our arguments for the S-wave results. The concluding 
section contains some discussion. 



2 Proof of the results 

The radial equation, for a two-body system with reduced mass [i in an 
attractive central potetial V(r) for u n ^{r) = rR n ^{r) is 

- Hu'U{r) + [Wi{r) - E n ]u ni (r) = 0, (5) 

where 
and 

The radial wavefunction R n ^{r) for energy E n is real so its modulus square 
is the same as its square. Consequently, the inequalities in the introduction 
are usually stated for the modulus square. The energy of the bound state 
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increases with principal quantum number, n, thus E\ < E2 < E3 . . . . The 
potential obeys the standard restrictions, namely, ]hn r ^,or V(r) = 0. Also, 
recall that R n i(r) behaves ELS V cLS r tends to zero. For an attractive force, 
the asymptotic behaviour (r — > 00) of the radial wavefunction u(r) will be 
like exp(— ar), a > 0. 

Multiply the radial equation by u' n p and integrate from zero to infinity. 
The term with E n gives zero. One integration by parts gives 

i>oo 

H[u' ne (0)}% = / Wl(r)u 2 ne (r)dr. (8) 
J 

The term We(r)u^(r) from the partial integration does not contribute. This 
is obvious for the upper limit r = 00. One has to be careful at the lower limit 
r = 0. However, since V(r) is less singular than r~ 2 and wtAr) ~ r 2(^+i) 
as r — > 0, the lower limit also does not contribute. All this is well known. 
Before specialising to S-wave it is interesting to consider the above equation 
for non-zero I. 

2.1 Result for non-zero I 

In this case, since the l.h.s. of Eq(8) is zero, the equation simply says that 
the expectation value of the effective force Wg(r) is zero. Alternatively, it 
implies that the expectation value of V'(r) for a general potential is related 
to that of r~ 3 . Explicitly, 

(V'(r)) ni = 2Hl(£ + l)(\) . (9) 

This general result (probably known personally to many [6]) deserves to be 
better known. It is is very useful. For example, for a Coulomb potential 
it immediately gives the correct relation between the expectation values of 
r~ 2 and r~ 3 . 

2.2 S-wave relations 

For £ = 0, Eq(8) reduces to 

/•oo 

HS n (0) = / V'(r)u 2 n0 (r)dr = (V'(r)) n0 . (10) 
J 

This is a well-known result and provides the basis for the arguments leading 
to the proof of the results given in the introduction. 
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Case (a): V"(r) = for all r. 

This is the case of the attractive linear potential Vr) = Xr. So, V'(r) = A 
is positive for all r. In this case, Eq.(lO) reduces to simply 

HS n (0) = A, (11) 

since u n o(r) is normalized, that is 

/•oo 

/ u 2 n0 (r)dr = l. (12) 
Jo 

Thus, in this case HS n (0) is a constant (the potential strength) independent 
of [i or n as required. It is well known that the linear potential is exactly 
solvable in terms of Airy functions. The above relation for HS n (0) has been 
noted earlier using the explicit solutions [7] . 

Case(b). V"{r) < and V'(r) > for all r, with V'(oo) finite. 

A well-known exactly solvable example of this case is the Coulomb po- 
tential. Perform an integration by parts in Eq(10) to obtain 

/•oo 

HS n (0) = y'(oo) - / V"(r)f n0 (r)dr, (13) 
J o 

where 




Note that / n o(oo) = 1 because the radial wavefunction is normalized . The 
term V'(r)/ n o(r), from the integration by parts at r = oo, gives V'(oo) 
while that at r = vanishes. This is because V(r) is less singular than 
r -2 as r tends to while one expects f n o( r ) ~ r 3 as r tends to because 
u no( r ) ~ 7-2 • Physically, / n o( r ) represents the probability of finding the 
particle (two-body system) between and r. As n increases the energy E n 
of the state increases. For a single well potential, the classical turning point 
E n = V(r) will be at a larger value of r as n or energy increases (keep in mind 
the Coulomb potential). So, one expects f n o( r ) to be larger than f m o(r) for 
n < m for small r. Asymptotically, the states with larger energy may have 
a longer tail (true for the Coulomb case) with the result / n o( r ) — / m o( r ) is 
negative for large r for n < m. However, the large r contribution will be 
damped out since for concave potential V"(r) decreases asymptotically. In 
fact, for power law potentials r k with 1 > k > —2 with V'(r) > 0, V"(r) 
goes to zero as r tends to infinity and goes to infinity as r tends to zero. 
Given the above physical reasoning, we expect that 

/•oo 

H[S n (0) - 5 m (0)] = - / V"(r)[f n0 (r) - f m0 {r)]dr > (15) 
Jo 
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for n < m if V"(r) < for all r. This gives the first inequality in Eq(3). 

For the variation with respect to the reduced mass, we note that with 
incresing /i, the system will shrink in size. So, physically one expects that 
fno(r) will increase, that is, [^[/no( r )] > (Sj. Thus, taking the derivative 
w.r.t. fj, of Eq(13 ), since V'(oo) is a constant, gives the second inequality 
in Eq(3). 

Case(c). V"{r) > and V'(r) > for all r, with V'(0) finite. 

An exactly solvable example for this case is harmonic oscillator potential 
r 2 . In this case we perform a slightly different integration by parts in Eq(10) 
to obtain 

POO 

HS n (0) = V\0) + I V"{r)g n0 {r)dr. (16) 

where 

" oo 

2 







9no(r) = / u n0 (r)dr = 1 - f n0 (r). (17) 



r 



The term V'(r)g n o(r), from the integration by parts, gives V'(0) for r = 0. 
For the upper limit r — > oo it vanishes because u n o(r) represents a bound 
state. From the above two equations we obtain 

/•oo 

H[S n (0) - 5 m (0)] = / V"(r)[g n0 (r) - g m0 (r)]dr (18) 



J o 

Now, as n or E n increases, g n o( r ) w iH increase so for a convex potential 
V"{r) > for all r, [S n (0) — S m (0)] > for n > m. For the variation with 
respect to the reduced mass /U, we obtain 

| (^.,(0)) < 0, (19) 

since the variation with \x of g n o( r ) is opposite in sign to that of f n o( r )- This 
concludes the proof of the S-wave results given in the introduction. 

A justification for the behaviour of fno{r) or g n o( r ) w.r.t. n required 
above, in Eq(15) and Eq(18), for our proof is provided by the Virial theorem. 
For S-wave levels,it states 

E n0 = {U(r)U, U(r) = V(r) + Kv'{r). (20) 

So, if U (oo) is finite,then an integration by parts gives [8] 

;>oo 

- [E n0 - E m0 ] = / U'(r)[f n0 (r) - f m0 {r)]dr > 0. (21) 
J o 
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For n < m, the l.h.s. is always positive, it implies that the integral is always 
positive. If U'(r) > for all r , one can physically understand the positivity 
of the integral because of the expected behavior of f n o( r ) discussed prior to 



For example, general power law potentials V(r) = — \X\r a with —2 < 
a < (which includes the Coulomb potential) satisfy all the conditons 
required for Case(b) and Eq(21).For these potentials the energy E n Q is al- 
ways negative. Furthermore, the behaviour of U'(r) ~ r a - 1 [ s milder than 
that of (— V"{r) ~ r a ~ 2 because of the extra derivative in the latter. In 
fact,[(— rV") — U'\ = ||Aa|r a_1 which is positive for all r.Thus,the positiv- 
ity of the integral in Eq(21) implies the posivity of the integal in Eq(15) for 
(-V"(r) > for all r. 

Guided by the above discussion, for Case(b) the already stated conditions 
on the potential, perhaps need to be supplemented with the conditions U (oo) 
finite and U'(r) > for all r. 

For potentials with U(0) finite, an integration by parts of Eq(20) will 
give 



for n < m. Thus, if U'{r) > for all r, the integral will be negative since 
g n o(r) increases with n as noted earlier. 

For example, power law potentials V(r) = \X\r a with a > 1, with all 
E n o > 0, satisfy all the conditions needed for Eq(18) and Eq(22). Per- 
haps, for Case(c) we should require the potential to satsfy additional cond- 
tions , namely U(0) finite with U'(r) > for all r. 

Even if the general condtions on the potential have to be supplemented(as 
per the above discussion) for the cases (b) and (c),the fact remains that the 
proof of the inequalities is for all n for a general class of potentials. 

3 Concluding remarks 

The proofs of the S-wave results presented above are not highly mathemat- 
ical but rely on the physical meaning of the quantities involved and how 
they are expected to change physically with energy of the bound state (or 
n) and the reduced mass ji. All known soluble examples of ttractive poten- 
tials with curvature of the same sign for all r support the results given in 
the introduction. Such potentials imply that the bound states lie in a single 
potential well. This is important for the physical arguments presented here. 
A potential with more than a single well cannot possbly have cuvature of 



Eq(15). 



- [E nQ - E m0 ] = - U'(r)[g n0 {r) - g m o(r)]dr > 0, 




(22) 
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the same sign everywhere. There are lot of solvable potentials for S-waves 
[9] without definite curvature for which the results given in the introduction 
may or may not hold. It would be an interesting challenge to find a counter 
example to the inequalities presented in the introduction. 
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